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Abstract — In this work the crystallisation of palm oil (a typical melt) is used as a model process and kinetic parameters at 
different end temperatures and cooling rates were evaluated using three model approaches: two from the traditional melt 
fractionation, one from the classic solution crystallisation. The objective was to establish critical understanding on the 
various model approaches and their applications. Differential Scanning Calorimetry (DSC) was used to obtain the melting 
temperatures, isothermal induction times and enthalpy of the crystallisation for the Fisher-Turnbull and Avrami model 
evaluations, while turbidity and temperature probes were utilized to generate metastable zone width as a function of cooling 
rates for the classical Nyvlt model analysis. Our results show that the Fisher-Turnbull and the Nyvlt models are useful in 
estimating the nucleation rate constants (kn) with reasonable agreement: this unites the model approaches and allows 
comparison between fat fractionation and solution crystallisation of organics. While the Avrami model is capable of 
evaluating the growth mechanism of the formed crystals (n) and the overall crystallisation rate constant (k), none of the 
parameters can be compared with other models because of the different definitions of the growth mechanism as well as the 
dependence of k value on the growth mechanism when dealing different fats of varying compositions. This is the first 
contribution of this work. In addition, our results indicate that the growth of nuclei to a stable size was generally slower in 
melt crystallisation compared to a typical organic solution system due to the relatively lower rate constant (kn) and the high 
viscosity and multicomponent properties of the melt system used in this study. This is the second contribution of this work. 
Keywords — Avrami model, Crystallisation, Fisher-Turnbull model, nucleation rate constant, Nyvlt model, palm oil 
fractionation. 

I. Introduction 

Crystallisation by cooling is a well-known unit operation used in a wide range of industries from pharmaceuticals to food for 
product purification, separation and/or recovery. For pharmaceutical ingredients, for example, the purity of the solute out of 
the solution is a critical product parameter. For edible oils/melts in food industry, on the other hand, it is the purity or the 
clearness of the solution that matters. Although the terminology of crystallisation is used for all inorganic and organic 
compounds, while fractionation is applied to all melts, the essence of the process is the same, i.e. by cooling for all. In spite 
of the similarities in the processes, different model approaches have been used for crystallisation from solution (e.g. APIs) 
than from melts (e.g. fats). Taking palm oil as the model compound for melt, it consists of a wide range of triglyceride 
(TAGs) fractions based on their physical state at different temperatures 111 - olein (liquid), mid-fractions (soft solid) and 
stearin (solid). These are significantly different to crystallisation of paracetamol for instance. The aims of this work are to 
evaluate and compare kinetics from a melt crystallisation extracted using models for both melts and organics; and to establish 
critical understanding on the various model approaches and their applications. 

1.1 The Models 

Three models are considered: the first two from traditional melt crystallisation, the last one from solution crystallisation. 

The first model for melts is the Avrami model dealing the rate of phase change (i.e. the formation of solid phase from liquid 
phase per unit time) at isothermal conditions. It is a simplistic model that uses the deterministic 121 approach to evaluate the 
overall kinetics in relation to nucleation mode and crystal growth mechanism 13 ' 41 of the isothermal process. The model was 
developed based on a few assumptions, one of which is that the growth rate is constant throughout the crystallisation 
process 131 ; but this may not be the case for fat systems in practice due to their multicomponent property. Hence modifications 
to the Avrami model have been proposed by other researchers [5 " 7] ; nevertheless the Avrami model, though simplistic, has 
widely and successfully been used to characterise the overall kinetics of fat crystallisation, 121 that is why it is chosen for this 
work. 

The second model in consideration in this study is the Fisher-Turnbull (F-T) correlation, which estimates the activation 
energy barrier required for the formation of stable nuclei by exploring the relationship between the degree of supercooling 
and the induction time needed for the formation of the nuclei 181 . Similar to the Avrami approach, there are a few assumptions 
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for the F-T model, e.g. induction time is a direct indication of nucleation event at a macroscopic level 19, 10] . This assumption 
is not unreasonable due to the negligibility of some of the inherent steps 1 " 1 . The activation energy would allow the nucleation 
rate constant at isothermal conditions to be determined. Although this method was originally designed for pine systems 
which is not the case for fat systems, it has widely and successfully been used to characterise the nucleation kinetics of fat 
crystallisation [12] ; hence it is adopted for this work. 

Since fractionation of fat is effectively a cooling crystallisation process, the classical Nyvlt nucleation theory is also 
applicable where the nucleation rate is a function of supersaturation and the rate of supersaturation generation is a function of 
cooling rate. By measuring the metastable zone width at different cooling rates experimentally, the apparent nucleation order 
and the nucleation rate constant can be determined and compared with these from the two holistic models. 

II. Materials and Method 


2.1 Materials 


Refined, bleached and deodorised (RBD) palm oil was supplied by AarhusKarlshamn, AAK (Hull, UK). The saturated fat 
content was 50.8% with palmitic acid forming 44% of the fatty acid assay. The same batch of material was used for all the 
experiments without further purification. Palm stearin and olein were considered as the two main fractions that constitute the 
solid and liquid products respectively. 

2.2 Isothermal Crystallisation 

The induction times of crystallisation (x) and the melting point (T,,,) were determined with a NETZSCH differential scanning 
calorimeter (model DSC 214 Polyma, NETZSCH, Germany) equipped with an intracooler cooling system operating between 
a temperature range of -70 to 600°C. Nitrogen at a flowrate of 40 mL/min and compressed air at a flowrate of 20 mL/min 
were used as the purge and protective gas respectively. 15-18 mg of sample was weighed into an aluminium pan and sealed 
with a pierced lid; an empty and sealed pan was used as a reference. The sample was heated to 353K at a rate of 10 K/min 
and held at this temperature for 10 min. The sample was then cooled to the pre-set end temperature (293 - 301K) at 20 K/min 
and held at each of the end temperatures for crystallisation to occur. Since the crystallisation is an exothermic process, the 
onset of the exothermic peak during the isothermal period was taken as the induction time (t) and all were evaluated using 
the DSC 214 software. At the end of the crystallisation process, the sample was then heated back to 353K at 10 K/min to 
obtain the melting curves. 

2.3 Non-isothermal Crystallisation 

To obtain the dynamic thermograms of the palm oil used for this study, the same DSC equipment as above was used; the 
sample (10 - 15 mg) was heated to 353K for 5 min, then cooled to 223K at the rate of 10 K/min; it was held at this 
temperature for 5 min to allow complete crystallisation. The sample was then heated again to 353K at the rate of 10 K/min to 
obtain the melting curve; the melting point (T m ) of the palm oil sample was taken as the temperature at the end of the curve, 
i.e. when melting is completed. This method of the determination of melting point is in accordance to previous works by 
Dodd & Tonge 1131 and Nassu & Goncalves" 41 . 

For the evaluation of the metastable zone width (A T max ) at different cooling rates, about 400 - 600 mL of palm oil was added 
to a crystalliser vessel - a DN 100 jacketed round-bottom stirred tank fitted with 5 -port glass lid and PTFE anchor impeller; 
the oil was melted at 50°C to ensure complete dissolution and then cooled to 20°C at different cooling rates (0.25 - 1.0 
°C/min) at constant mixing condition. GP200-R2 water bath (Grant Inst. Cambridge Ltd) was used to control the temperature 
of the vessel jacket. Turbidity and temperature (PT100) probes (CrystalEyes system from HEL, UK) were inserted into the 
crystalliser vessel to monitor and record the cloudiness and temperature of the solution. The system was interfaced with a PC 
for real-time display, logging and data analysis. 

2.4 Estimation of kinetic parameters 

The Avrami model. The basic principle can be illustrated by imagining raindrops falling in a puddle. The raindrops produce 
expanding circles of waves which intersect and cover the whole surface. The drops may fall sporadically or all at once. In 
either case, they must strike the puddle surface at random points. The expanding circles of waves are the growth front of the 
spherulites, and the points of impact are the crystallite nuclei. The original derivations by Avrami [3] have been simplified by 
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Evans 1151 and put into polymer context by Meares 1161 and Hay 1171 . Through probability derivations 1181 , the volume fraction of 
crystalline material, X, known widely as the degree of crystallinity, can be written as: 

l-X =e~ E (la) 

Where E represents the average number of fronts of all such points in the system. For low degrees of crystallinity, a useful 
approximation is X » E. For the bulk crystallisation of polymers, X in the exponent may be considered the volume or volume 
fraction of crystalline materials, V t , i.e.: 

l-X=e~ V ' (lb) 

This has been the widely accepted and used equation in bulk crystallisation 119 ’ 221 . For either instantaneous or sporadic 
nucleation, equation (lb) can be written as: 

l — X=e Kt or — ln(l — x) = kt n (lc) 

where x is the fraction of solid formed at a particular time (t)\ k the overall crystallisation rate constant; and n the index of 
crystallisation, also referred to as the Avrami exponent, the phenomenological index of crystallisation 1231 , depending not only 
on the structure of the crystal, but also on the nature of nucleation 14 ' 24] . For example, when n = 1, it corresponds to rod-like 
growth from instantaneous nuclei; whereas n = 3 or 4 refers to spherulitic growth from either sporadic or instantaneous 
nucleation 1 " 81 . However, fractional values of n also exist due to secondary crystallisation, e.g. lower n values (<1) are caused 
by linear crystal growth 1261 . 

The value of x was estimated by integrating the isothermal crystallisation exothermic peak from the DSC measurement, as 
given by the eq. 2 below where AH, is the heat of crystallisation (W/g) which is evaluated as the area under the crystallisation 
curve from the onset of the peak to time t using the DSC 214 software; All,,,, is the total area under the curve. More detailed 
description of the methodology can be found elsewhere 1271 . The kinetic parameters, n and k, were obtained from the slope and 
the intercept of the plot of [ln(-ln(l-x)] against ln(f) respectively. The unit of k here is the inverse of time to the power of n, 
depending on the index of crystallisation: 



The Fisher-Turnbull correlation: This model which was originated by Becker[28] for condensed systems (liquid-solid or 
solid-solid transformations) links nucleation rate with the activation free energy required for the formation of a stable nucleus 
and the diffusion energy for phase change. The former (AGc) was estimated in the equations below: 
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Where h and k b are the Planck (6.626xl0 34 Js) and the Boltzmann (1.381 x 10" 23 JK" 1 ) constants respectively; N and N A are 
the Avogadro’s number and constant (6.022 x 10"’ and 6.022 x 10 23 mol" 1 respectively); AH is the enthalpy of fusion (Jmol" 
' ); T m the melting point of the palm oil (K); AT the degree of supercooling (T m - T) with T being the crystallisation 
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temperature; s the slope of the plot of In (Tx) against 1/T (AT) 2 with the unit of K 3 ; AGd the activation energy for molecular 
diffusion. Using the activation energy (AG C ) from this model, the nucleation rate constant (k n ) can be estimated based on the 
Arrhenius -type equation 1291 shown below, where R is the molecular gas constant (8.314 Jmol^K" 1 ): 


k 


n 



(4a) 


If the rate of solid crystals formation ( dx/dt ) is an indication of the crystallisation rate (J). then the rate equation is given as: 


J 



(4b) 


Assume 


f(.x) = (l - x) 


(4c) 


where fix) is referred as a function dependent on the kinetic model of the crystallisation process 1301 . If homogeneous 
nucleation is assumed, the f(x) is taken as (1 - x)" where n is the nucleation order 1291 . Hence using the DSC data obtained, a 
plot of ln(dv/dr) against ln(l-x) will give a slope equal to n, and the value of A can be obtained from the intercept and the AG, 
evaluated from eq. (3c) for the different end temperatures considered. Using these A values, the rate constant (k„) can 
therefore be estimated according to eq. (4a). The unit of k n is the inverse of time (s’ 1 ). 

The classic Nyvlt model. This model is based on the assumption that when nuclei is detectable, the rate of supersaturation is 
equal to the rate of nucleation. From the classical nucleation theory, the primary number nucleation rate, J n , can be expressed 
as a function of supersaturation, AC 1111 : 

J n = k n (AC nlax ) (5a) 


Where k n is the nucleation rate constant on a number basis, AC the supersaturation and n the nucleation order. The unit of J„ 
is the number of crystals generated per volume and per time. For cooling crystallization, the rate of supersaturation 
generation is a function of the cooling rate, /?: 


Jn=fi 

V 


dC sat 

clT 


\ 


J 


(5b) 


Where dC sat /dT is the slope of the solubility curve for a given saturation temperature. For the palm oil used in this work, the 
solubility curve of stearin/olein as a function of temperature ( T in °C) is given by C sat = lxlO" 6 e 0,2695T . At the point of 
nucleation, the maximum possible supercooling (or metastable zone width (MSZW)), AT max , and the corresponding 
maximum supersaturation, AC max , are given by: 
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Combining eq. (5b) and (5c) into eq. 5a results in: 
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The Nyvlt model effectively relates the effect of cooling rate (/7) with the metastable zone width ( AT max ) at a particular 
operating condition as (see Nyvlt 1311 and Sangwal 1321 for full derivation of this model): 
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As a result from a plot of ln(j3) against ln(AT max ), the nucleation order («) and the nucleation rate constant (k n ) can be 
estimated from the slope and the intercept of the plot respectively. For number basis, the unit of k n is the inverse of time (s' 1 ). 

In summary, the Avrami model evaluates the overall crystallisation rate constant (s" n ) and the crystallisation growth 
mechanism at isothermal conditions; the F-T model estimates the nucleation rate (s' 1 ) from the activation energy required for 
nucleation at a given temperature and finally, the Nyvlt model predicts both the nucleation rate (s' 1 ) and the nucleation order 
at non-isothermal conditions. 


III. Data for models 

For the palm oil, the solid phase is the stearin while the liquid phase is olein, which are represented in the cooling and 
melting thermograms (DSC output - the indication of enthalpy change as a function of temperature) obtained from the 
dynamic DSC analysis. The cooling thermogram (Fig 1) shows two distinct exothermic peaks: the one at a low temperature 
(-4.5°C) corresponds to the low melting fraction (olein) due to its high degree of unsaturation; while the second peak at 16°C 
is associated with the crystallisation of the high melting fractions (stearin). The broader peak at this low temperature than at 
the high temperature indicates a relatively wider range of triglycerides within the temperature region. Similar results were 
reported in a previous study on RBD palm oil 1331 . 

Likewise, the melting thermogram (Fig 2) also shows two groups of endothermic peaks; the peaks at 28°C and 35.5°C 
correspond to the high melting fractions (stearin) while the low temperature endothermic peaks (4.5°C and 7.5°C) correspond 
to the low melting fractions (olein); the larger size of this endotherm is consistent with the wider range of low melting 
fractions discussed earlier. 



Fig 1: Dynamic cooling thermogram ol palm 

OIL AND ITS TRACTIONS 



Temp (°C) 


Fig 2: Dynamic melting thermogram ol palm 
oil and its tractions 


The isothermal crystallisation thermograms at different end temperatures from 293 K to 301 K are shown in Fig 3 where the 
heat of crystallisation is plotted as a function of time. It can be seen that the time taken (i.e. induction time, r) for 
crystallisation to begin - taken as the onset of the exothermic peak - increased with increasing end temperatures. This can be 
attributed to reduced degree of supercooling as the higher the end temperature, the closer to the equilibrium temperature the 
system becomes, hence the lower the driving force for the crystallisation process, as it is case for solution crystallising 1141 . 
Images from the polarised light microscope in Fig 4 show different crystal morphologies from the samples collected after 
about 1 hr for the different crystallisation temperatures; for instance, fully-formed spherulitic crystals were observed at 293 
K, dendritic spherulites at 303 K and a developing spherulites at 298 K, which supports only one morphology at a given 
crystallisation temperature studied. 
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Fig 3. Isothermal crystallisation thermogram showing the induction time at different end 

TEMPERATURES 



Fig 4: Microscopic images of palm stearin crystals at 293 K (a), 298 K (b) and 303 K (c) 

IV. Discussion 


4.1 The Avrami model 

Using the DSC isothermal crystallisation thermograms and the induction time analysis, a plot of the solid fraction (x) against 
time it) is shown in Fig 5 with the well-known sigmoid curve for isothermal fat crystallisation. The changing gradient of the 
curve depicts the initial formation of nuclei, followed by a rapid crystal growth, then the slowing down by impingement due 
to colliding crystal faces [2 ’ 3SI , as indicated on the right of Fig. 5. These stages are relatively more defined as “5” shaped 
curves at the higher temperatures (>296 K) due to the faster crystallisation rate at low temperatures in a previous study by 
Zhang et a/ 1361 . 



Fig 5: Percentage of solid fat fraction with time at different crystallisation temperatures 
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By plotting ln[-ln(l-;t)] against ln(f) from the Avrami analysis, x values between 25% and 75% [41 are displaced as a function 
of the crystallisation time in Fig 6 a good linear relationship ( R 2 > 0.99) indicating the applicability of the model to the fat 
system used in this study. The resultant n and k values are listed in Table 1. It can be seen from this table that the Avrami 
exponent ( n ) for the palm oil used in this study was relatively unaffected by the end crystallisation temperatures within the 
range considered, as it was relatively constant (2. 5+0.2, P > 0.05 using ANOVA with single factor analysis). The n value is 
an indication of the growth mechanism of the crystals formed according to Kawamura [37] and Wright et a/ [3S| , which ranges 
from rod-like for n = 1 to polyhedral growth (three-dimensional) for n = 4. For polymer and organics systems, an n value of 2 
or 3 indicates one or two dimensional nucleation of the crystal nucleus; a range of 3.85 to 4.49 has been reported for 
polyethylene oxide. From n = 3 in this study, the crystal growth mechanism can be taken as plate-like (i.e. two-dimensional 
growth). However previous studies have reported n values greater than 4 for palm oil crystallisation which was attributed to 
sporadic and heterogeneous nucleation. This discrepancy could be caused by different triglyceride compositions and different 
cooling rates of the DSC analyses used (unpublished work by Toro-Vazquez et a/ [39] ). For the estimated overall rate constant 
( k ) also shown in Table 1, it generally decreases with increasing end temperatures, indicating a slower crystallisation process 
at higher temperatures as expected 1341 . The sharp increase in the k value from 296 K to lower temperatures can be linked to 
polymorphic transition occurring at those temperatures. In summary, the rate constants (k) obtained from this work fall 
within the similar range as the previous work, suggesting similar crystallisation kinetics . 


4 - 0.5 


- 1.5 

5 6 7 8 9 

ln(t) 

Fig 6: Plot for Avrami model analysis 

Note that the ^-dependency of the unit of these k values (s' n ) makes it impossible to compare with related values from other 
models; hence the trend of k with respect to the end temperatures, as opposed to the actual values, is the key element from 
this model. 

4.2 The Fisher-Turnbull model 

The results discussed above give an indication of the overall crystallisation kinetics. In order to examine the nucleation 
occurrence, the Fisher-Turnbull equation was used to estimate the activation free energy barrier (AG C ) needed for the 
formation of stable nuclei. The related plot of (ln(7r) against 1/T(AT) 2 ) and the calculated values of AG C together with the 
estimated nucleation rate constants ( k n ) using eq. (4a) are shown in Fig 7 and Table 1 respectively. Note that the melting 
temperature ( T,„ ) used to calculate AT was taken as the peak end temperature of the melting curve shown in Fig 2 which was 
44°C. Fig 7 effectively displays the relationship between the energy barrier for nucleation (the vertical axis) and the inverse 
of supersaturation (the horizontal axis); the higher the crystallisation temperature, the larger the energy barrier (higher AG C ), 
again due to lower degree of supercooling at higher end temperatures, which is as expected 1411 . 

Fig 7 has a good but discontinued linearity from about 298 K. This discontinuity can be attributed to the formation of crystals 
of different polymorphic states; according to Ostwald’s rule of stages 1421 , nucleation of the less stable polymorph (a) was 
favoured especially at low end temperatures (higher supersaturation) where nucleation rate is higher with less energy 
barrier 1431 ; hence it can be inferred that the crystals formed at T > 298 K is the P polymorph 1441 ; this trend is similar to results 
shown in previous studies by Chen et a/ 1451 and Ng 1461 . The features of the discontinuity and the different nucleation 
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mechanisms leading to different polymorphs are also clearly shown in Table 1 where higher nucleation rate constants (k„) 
were extracted at lower end temperatures and vice versa; the corresponding n values (see eq. 4c) obtained here was ~1 hence 
the nucleation process can be taken as a first order. 


Table 1 

Kinetics parameters for all Models; (Tm = 317K) 



Avrami model 

Fisher-Turnbull model 

Nyvlt model 

T C (K) 

n 

kx 10 10 (s' n ) 

A G c (kJ/mol) 

k„ x 10 6 (s' 1 ) 

n 

k„ x I0 6 (s~') 

293 

2.7 + 0.2 

2866 

1.37 

4188 

8.97 

399 

295 

2.3 +0.1 

1646 

1.63 

4158 

296 

2.2 + 0.0 

282 

1.79 

669 

297 

2.5 +0.1 

12 

1.98 

465 

299 

2.3 +0.4 

12 

4.75 

476 

301 

2.7 + 0.4 

4 

6.01 

202 


The property, AGd, being the activation energy for molecular diffusion takes into account the impeded molecular movement 
due to increased viscosity at low temperatures; hence it is expected that the nucleation rate becomes more influenced by AG d 
as the end temperature is reduced 1471 . This was however not explored further in this study. 



l /( AT) 2 Tx 10 6 


Fig 7: Kinetics parameter from Fisher-Turnbull equation 
4.3 The Nyvlt model 

To extract non-isothermal kinetics from the Nyvlt model, the measurements of metastable zone width (MSZW) as a function 
of cooing rates are required. Fig 8 shows the typical turbidimetry measurement, with the highest transmittances indicating 
that the palm oil is completely melted. The moment when the transmittance starts to drop is regarded as the onset of 
nucleation (i.e. the formation of palm stearin solid crystals). The transmittance keeps decreasing as crystallisation proceeds 
with more stearin crystals in the system until the signals level off to zero or close to zero indicating the end of the 
crystallisation process. The ‘dissolution’ profile was obtained by reheating the system from the end temperature (20°C) back 
to the start temperature (50°C). In this way, the metastable zone width (MSZW) can be determined as the difference between 
the ‘dissolution’ and the crystallisation temperatures at each of the cooling rates. Table 2 lists the MSZW for the three 
cooling rates investigated; as expected that the highest cooling rate (1.0°C/min) has the lowest nucleation temperature hence 
the largest MSZW due to rapid supercooling and longer nuclei stability period 148 ' 5#] ; similar results were reported in previous 
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studies relating to the effects of cooling rate on fat crystallisation 11, l This is also consistent with the results from the DSC 
isothermal analysis (Fig 3) discussed earlier in that the crystallisation rate is higher at lower end temperatures. 

By plotting ln(j3) against ln(A T max ), the nucleation order («) and the nucleation rate constant (k n ) can be evaluated according 
to eq. (5e) and are also shown in Table 1. Values of n between 6.0 and 16 were reported for fat systems of different 
tripalmiltin concentrations in a previous study 13 , and the value derived from this study (~9) is within this range. The k n value 
obtained here (3.99 x 10" 4 s" 1 ) is a lot lower than those reported for solution crystallisation of organics (e.g. adipic acid) in a 
previous study 1 ' 21 , and this may be attributed to the relatively slower rate of molecular growth to a stable nuclei for 
triglyceride systems due to the high viscosity; a higher nucleation order (~9) also supports the influence of mass transfer on 
melt crystallisation 1531 . 


Table 2 

MSZW VALUES FOR DIFFERENT COOLING RATES 


P (°C/min) 

T cry (°C) 

T dis (°C) 

AT max (°C) 

1.0 

24.2 + 0.2 

45.1+0.1 

21.0 + 0.2 

0.5 

25.2 + 0.3 

19.9 + 0.4 

0.25 

27.2 + 0.4 

17.9 + 0.4 


OJ 

O 


E 

CO 



Fig 8. Turbidity measurements for MSZW determination at different cooling rates 


4.4 Comparison of models 


Comparing the nucleation rate constants evaluated between the Fisher-Turnbull and the Nyvlt models (see in Table 1), the 
Nyvlt non-isothermal nucleation rate constant is of similar order of magnitude to the Fisher-Turnbull isothermal ones for the 
higher end temperatures (298 - 300 K), while of an order of magnitude lower for the lower end temperatures. The values 
themselves match reasonably well for the cooling rates investigated, giving the fact that the model assumptions are widely 
different, where the Nyvlt model took a polythermal approach to estimating the kinetic parameters, assuming that 
supersaturation rate corresponds to nucleation rate from the onset of nuclei formation and there is no growth in nuclei [s4] , 
while the Fisher-Turnbull model was based on isothermal conditions with homogenous nucleation being assumed in the 
crystallisation process. The two models can be used for estimating nucleation parameters for both melts and organics. 

The dependency of the k values on the growth mechanism (n) evaluated from the Avrami model together with the fact that it 
describes the overall crystallisation rate makes any comparison irrelevant. Furthermore the growth mechanism (n) from the 
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Avrami model has a completely different definition from the apparent nucleation order predicted from the Nyvlt model 
which is dependent on the formation process of stable nuclei 1321 ’ once again no comparison is possible. 

V. Conclusion 

In this work we use the crystallisation of palm oil (a typical melt) as the model process and evaluated the kinetic parameters 
using three model approaches: two from the traditional melt fractionation, one from the classic solution crystallisation. The 
objective was to establish critical understanding on the various model approaches and their applications. Our results show 
that the Fisher-Turnbull and the Nyvlt models are useful in estimating the nucleation rate constants (k„) with reasonable 
agreement: this unites the model approaches and allows comparison between fat fractionation and solution crystallisation. 
The lower rate constants predicted from the melt crystallisation in comparison to that from the solution crystallisation of 
organics are expected and resulted from the high viscosity that impedes the growth of nuclei to a stable size. 

While the Avrami model is capable of evaluating the growth mechanism of the formed crystals (n) and the overall 
crystallisation rate constant ( k ), none of the parameters can be compared with other models because of the different definition 
of the growth mechanism as well as the dependence of k value on the growth mechanism when dealing different fats of 
varying compositions. 

Finally, the end temperatures and cooling profiles are the essential operational parameters that can be explored in order to 
obtain crystal product of desired physical properties. 
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